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Abstract 

In this paper we describe an algorithm for implicitizing rational hyper- 
surfaces in case there exists at m ost a finite number o f base points. It is 
based on a technique exposed in Buse and Jouanoloul |2002| |. where im- 
plicit equations are obtained as determinants of certain graded parts of 
a so-called approximation complex. We detail and improve this method 
by providing an in-depth study of the cohomology of such a complex. In 
both particular cases of interest of curve and surface implicitization we 
also yield explicit algorithms which only involves linear algebra routines. 



1. Introduction 

The implicitization problem asks for an implicit equation of a rational hyper- 
surface given by a parameterization map <fi : P n_1 — > P n ; with n > 2. This 
problem received recently a particular interest, especially in the cases n = 2 
and n = 3, because it is a key-point in computer aided geometric design and 
modeling. There are basically three kinds of methods to compute such an im- 
plicit equation of rational curves and surfaces. The first methods use Grobner 
bases computations. Even if they always work, they are known to be quite slow 



in pra ctice and are hence rarely used in geometric modeling (see e.g. iHoffmann 
|l989l |). The second methods are based on resultant matrices. Such methods 
have the advantage to yield square matrices whose determinant is an implicit 
equation. This more compact formulation of an implicit equation is very useful 
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and well behaved algorithms are known to work with . Howe ver such methods 



sec 



Jouanolou 19961), and in case base 



are only known if there is no base point 
points are isolated and local co mplete intersection with some additional techni- 
cal hypothesis (see iBusel |200l| ) . The third meth ods are based on the syzyg ies 



as 



of the parameterization. They were introduced in iSederberg and Chen 
the method of "moving surfaces". For curve implicitization this method, also 
called "moving lines", reveals very efficient and general. This generality is no 
longer tru e for s urfac e implicitization, but the method remains very efficient. In 



Cox et al.l |2000| and ID 'Andrea! |200I| its validity was proved in the absence of 



base p oint. An extension in the presence of base points is exposed in lBuse et al 



[2002] , assuming five base points hypothesis, one being that the base points are 
isolated and locally complete intersection. 



A s imilar approach involving syzygies was recently explored in lBuse and Jouanolou 
20021 ] using more systematic tools from algebraic geometry and commutative 
algebr a, among them the so-called approximation complexes (see IVasconcelosI 
" l994j ). A new method was given to solve the hypersurface (and hence curve 



and surface) implicitization problem in case base points are isolated and locally 
complete intersection, without any additional hypothesis. 



In t his paper we detail and improve the results exposed in lBuse and Jouanolou 
|2002j | for implicitizing rational hypersurfaces in case the ideal of base points is 
finite and a projective local almost complete intersection. We also give some 
other properties based on the exactness on these so-called approximation com- 
plexes. Section 2 and section 3 provide a complete and detailed description of 
our results in the particular cases of interest of respectively curve and surface 
implicitizations. Algorithms are completely describe and we present some exam- 
ples. Section 4 deals with the general case of hypersurface implicitization, and 
is much more technical that both previous one. In particular all the proofs are 
given in this section. 

Hereafter IK denotes any field. 



2. Implicitization of rational parametric curves 

In this section we present a method to implicitize any parameterized plane curve. 
Let /c/1,/2 be three homogeneous polynomials in K[s, £] of the same degree 
d > 1, and consider the rational map 

(s:t) h- (f (s,t) : fx(s,t) : f 2 (s,t)). 

We denote by x,y,z the homogeneous coordinates of P 2 . The algebraic closure 
of the image of is a curve in P 2 if and only if <ft is generically finite onto its 
image, that is 5 := gcd(/o, fi, $2) is not of degree d. We assume this minimal 
hypothesis and denote by C the scheme-theoretic closed image of 0. It is known 
that C is an irreducible and reduced curve of degree (d — deg(<5))//3, where (3 
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denotes the degree of onto its image (that is the number of points in a general 
fiber). Thus any implicit equation of C is of the form P(x,y,z) where P is an 
irreducible homogeneous polynomial of degree (d — 5)/ (5. 

In what follows we describe an algorithm to compute explicitly an implicit 
equation P of C (in fact we will compute P@) without any other hypothesis on 
the polynomials f , A, f 2 . I n the case deg(5) = we recover the well- known 



method of moving; lines (see ISederberg and Chenl |l99^ . ICox ei~aT1 |l99«| ). 



2.1. The method 

Let us denote by A the polynomial ring K[s,t]. We consider it with its natural 
graduation obtained by setting deg(s) = deg(t) = 1. From the polynomials 
fo, A> A of the parameterization 0, we can build the following well-known graded 
Koszul complex (the notation [— ] stands for the degree shift in A) 

-> A[-3d] ^ A[-2df % A\-df ±* A, (1) 

where the differentials are given by 

/ A \ (fx A \ 

d 3 = -A ],d 2 = -/„ / 2 ,di = (/ A A)- 

\fo J \ o -A -A/ 

In what follows we will not consider exactly this complex, but the complex 
obtained by tensorizing it by A[x,y,z] over A. This complex, that we denote 
{K,(f , A, A), «•)» is of the form 

- A[x, y, z] [-3d] ^ A[x, y, z] [-2d] 3 ^ A[x, y, z] [-df ^ A[x, y, z], 

where the matrices of the differentials di and ui are the same, for alH = 1, 2, 3. 
Note that the ring A[x, y, z] is naturally bi-graded, having a graduation coming 
from A = K[s, t], and another one coming from K[x, y, z] with deg(x) = deg(y) = 
deg(z) = 1; we hereafter adopt the notation (— ) for the degree shift in K[x, y, z\. 

We form another bi-graded Koszul complex on A[x, y,z], the one associated 
to the sequence (x,y,z). We denote it by (K,(x, y, z), v m ), it is of the form 

- A[x, y, z] (-3) ^ A[x, y, z] (-2) 3 ^ A[x, y, z] (-1) 3 ^ A[x, y, z], 

and the matrices of its differentials are obtained from the matrices of the dif- 
ferentials of fllj by replacing f by x, A by y and f 2 by z. Observe that since 
(x,y,z) is a regular sequence in A[x,y,z], the previous complex K,(x,y,z) is 
acyclic, that is to say all its homology groups Hi(K,(x, y, z)) vanish for i > 0. 

We can now construct a new bi-graded complex of A[x, y, z]-modules, denoted 
Z,, from both Koszul complexes (K m (f , A, A)> u ») an d (K»( x , z )i v ») (observe 
that these complexes differ only by their differentials). Define := ker(dj) for all 
i = 0, . . . , 3 (with do : A — >• 0), and set Zi := Zi[id] ®a A[x, y,z] for % = 0, . . . , 3, 
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which are bi-graded A[x, y, z] -modules. The map v\ induced obviously the bi- 
graded map, that we denote also by v±, 

Zi(-l) ^ Z = A[x, y, z] 

[91,92,93) >-> gix + g 2 y + g 3 z. 

Using the differential v 2 of K,(x,y,z) we can map Z 2 to A[x, y, z](— d) 3 , but 
since u\ o v 2 + v\ o U2 = (which follows from a straightforward computation), 
we have v 2 (Z 2 ) C And in the same way we can map Z% to Z 2 with the 
differential V3, since «2 o t> 3 + 1> 2 o w 3 = 0. Thus we obtain the following bi-graded 
complex (it is a complex since (K,(x, y, z), v m ) is) : 

(Z.,v.) : - Z 3 (-3) ^ Z 2 {-2) ^ Z 1 (-l) ^ Z = A[x,y,z]. 

This complex is known as the approximation complex of cycl es associated to the 
)olvn omials fo, fx, f 2 in K[s, t]. It was originally introduced in lSimis and Vasconcelosl 
1981] for studying Rees algebras through symmetric algebras (see also I Vasconcelosl 



1994]). 



Remark: Point out that in the language of the moving lines method, the image 
of a triple (gx, g 2 , #3) G ^iM(o) by V\ is nothing but a moving line of degr ee v 



follow ing the surface parameterized by the polynomials fo, fx, f2 ( see e -g- \Cox 
|20m| ). 

Since we have supposed that f , fx, f 2 are all non-zero we have Z 3 = 0, and 
consequently Z 3 = 0, and the following theorem (recall 5 := gcd(/o, fx, ^2)) : 

Theorem 2.1: The determinant of the graded complex of free K[x, y, z]-modules 

-> 2: 2[d _i](-2) ^ £![,,_!] (-1) ^ £ [d-i] = 

zs P(x,y, z)P , where P is an implicit equation of the curve C. 

Moreover, if deg(<5) = then Z 2 u_x] — 0; thus the dx d determinant of the 
map Z 1[d „ i](-l) ^* 2b[d-i] egwa/s P(x,y,z) . 



Proof: See section 5 in Buse and Jouanoloul 2002^ for a proof. We only mention 
how the last statement follows from the description of the approximation complex 
Z,. If deg(<5) = then depth( s ^(/o, fi, ^2) = 2, which means that fo, fx, f 2 have 
no base points in P 1 . This implies that not only the third homology group of the 
Koszul complex (JTJ vanishes, but also the second. It follows that Z 2 ~ A[— 3d], 
and hence Z 2 ~ A[x, y, z] [—d]. □ 

The second statement of this theorem gives exactly the matrix constructed 
by the method of moving lines. The first statement show that this method can 
be extended even if we do not assume deg(5) = 0, P 13 b eing obtained as the 
quotie nt of two determinants of respective size d and 5. In Buse and Jouanoloul 
|2002j | it is in fact proved that for any integer u > d — 1 the determinant of the 
complex (Z m )[ v ] equals In case deg(<5) = 0, this and the graded isomorphism 
Z 2 ~ A[x,y, z][— d] explain clearly, in our point of view, why the method of 
moving lines works so well with and only with moving lines of degree d — 1 . 
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2.2. The algorithm 

We here convert theorem 12.11 into an explicit algorithm where each step reduces 
to well-known and efficient linear algebra routines. 

Algorithm (for implicitizing rational parametric curves): 

Input : Three homogeneous polynomials f (s, t), fi(s, t), f 2 (s, t) of the same 

degree d > 1 such that 5 = gcd(/ , fi, $2) is n °t of degree d. 

Output : Either : 

• a square matrix Ai such that det(Ai) equals P 13 , in case deg(<5) = 0, 

• two square matrices Ai and A 2 , respectively of size d and deg(5), such that 
£^ equals 

1. Compute the matrix F-i of the hist map of (JIJ : Aj_ t ^2d-i- Its entries 
are either or a coefficient of /o, fi or / 2 , and it is of size 2d x 3d. 

2. Compute a kernel matrix Ki of the transpose of Fi . It has 3d columns and 
rank(^iw_ 1 i) lines. 

3. Construct the matrix Z t by Z 1 (i, j) = xKi (j, i) + yKt (j, i + d), zKj, (j, i + 2d) , 
with i = 1, . . . , d and j = 1, . . . , rank(^ 1 [ d _ 1 j). It is a matrix of the map 

Zi[d-i\{— 1) 2o[d-i]- 

4. If Zi is square then set Ai := Z^ else 

(a) Compute a list L\ of d integers indexing d independent columns in Z 1 . 
Let Ai be the dx d submatrix ofZ^ obtained by removing columns not 
in L\. 

(b) Compute the matrix F 2 of the second map of (JTJ) : A\_ x -Ajjd-iJ an< ^ 
a kernel matrix K 2 of its transpose. The matrix K 2 has 3d columns and 
rank(Z 2 [ d _!]) lines. 

(c) Construct the matrix Z' 2 by, for all j — 1, ... , rank(^ 2 w-i]); 

i = l,...,d : Z' 2 (i,j) =yK 2 (j,i) + zK 2 (j,i + d), 

Aihj) = ~xK 2 (j, i-d) + zK 2 (j, i + d), 
z 2(hj) = -xK 2 (j,i- d) -yK 2 (j,i). 



i — d + 1, . . . , 2d 
i = 2d + l,...,3d 



(d) Construct the rank(Z 1 [ d _ 1 j) x rank(^ 2 [rf-ii) matrix Z 2 whose j th column 
^(•jj) is the solution of the linear system t Z 2 (»,j).K 1 = Z 2 («, j). It is 
a matrix of the map Z2[d-i](— 1) Zi[d-x\- 

(e) Define A 2 to be the square submatrix ofZ 2 obtained by removing the 
lines indexed by L\ . 

endif 

Remark: In order to keep this algorithm easily understandable we kept the sym- 
bolic variables t and z, but they are of course useless and should be specialized 
to 1. 
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2.3. An example 

As we have already said, the previous algorithm is exactly the well-known method 
of moving lines in case deg(S) = 0. The only thing new is its extension to the 
case where deg(5) > 0. We illustrate it with the following (very simple) example. 

Let fo(s,t) = s 2 , fi(s,t) = st and f2(s,t) = t 2 . Applying our algorithm we 
found that the matrix Zi (step 3) is square : 

-y ~z 

x y 

We deduce that an implicit equation is given by xz — y 2 . Now let us multiply 
artificially each /j by s. We obtain new entries for our algorithm which are : 
fo(s,t) = s 3 , fi(s,t) = s 2 t and f?.{s,t) = st 2 . In this case the matrix Z t is not 
square, it is the 3x4 matrix : 

' —y —z —z 
Zi = | x y —z 
x y 

f-y -z -z x 

The matrix Ai can be chosen to be the submatrix of Z t given by x y 

\ x 

with determinant — x 2 z + xy 2 . Continuing the algorithm we obtain the matrices 



Z = 



y 
o 



, and A 2 = (— x). It follows that an implicit equation is — xz + y 2 . 



\-xJ 



3. Implicitization of rational parametric surfaces 

We arrive to the much more intricate problem of surface implicitization. Let 
foifiifiifz De f° ur homogeneous polynomials in K[s, t, tt] of the same degree 
d > 1, and consider the rational map 

(s : t : u) ^ (/ (s, t, u) : t, u) : / 2 (s, t, u) : / 3 (s, t, u)). 

We denote by x,y,z,w the homogeneous coordinates of P 2 . The closure of the 
image of is a surface if and only if the map is generically finite onto its 
image, which we assume hereafter. Thus let S denote the surface in P 3 obtained 
as the closed image of 0. Our problem is to compute an implicit equation of 
S. It is more difficult than curve implicitization because of the presence of base 
points in codimension 2 arising from the parameterization, base points in codi- 
mension 1 being easily removed by substituting each /j, % — 0, . . . , 3, by fi/5, 
where 5 := gcd(/oj fi, f2, fs), if necessary. Consequently, we assume that the 
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ideal I = (/o, /i, /2, /a) in K[s, £, tt] is at least of codimension 2, that is defines 
n isolated points j»i , . . . , p„ . in P 2 , and we have (see section 2 and theorem 2.5 in 
Buse and Jouanoloul 2002 ] for a proof and more details) : 



THEOREM 3.1: Let e Pi denote the algebraic multiplicity of the point for all 
% = 1, . . . ,n, and (3 denote the degree of <f> onto its image. Then 

^2 y- f /3deg(«S) if 4> is generically finite, 

^ Pi \ if (J) is not generically finite. 

i=X v - 

This theorem gives us the degree of the closed image of <fr if it is a surface (that 
we have denoted by S), and also that this image is not a surface if and only if 
=i e Pi- We will call an implicit equation ofS an equation of its associated 
divisor, which is an irreducible and homogeneous polynomial P(x, y, z, w) of total 
degree (d 2 - EILiO// 3 - 

In the following subsections we mainly present a method based on the approxi- 
mation complexes to compute explicitly an implicit equation of S (in fact we will 
compute P 13 ) in case the ideal / is a local complete intersection of codimension at 
least 2, and we also give some other related results. Then we expose an explicit 
description of the algorithm it yields and illustrate it with some examples. 



3.1. The method 

We denote by A the polynomial ring K[s, t, u] which is naturally graded by 
deg(s) = deg(t) = deg(w) = 1. Let us form the Koszul complex of fo, fi, f2, f3 
in A : 

-> A[-Ad] ^ A[-3d] 4 ^ A[-2df ^ A[-of ^ A, (2) 
where the differentials are given by 



f-h\ 

h 

-fx 
\foJ 



( h 


h 





o \ 


-h 





h 





fo 








h 





-fx 


-h 








fo 





-h 


V o 





fo 


h) 



( 


-fx 


-h 





-h 





\ 




fo 





-h 





-/a 










fo 


fx 








-h 


\ 











fo 


fx 


h) 



di=(/o fx h h). 



As for curve implicitization, we denote by (K,(f , fx, f 2 , fz),u.) this Koszul com- 
plex tensorized by A[x\ := A[x, y, z, w] over A, which is of the form : 



-> A[x][-4d} ^ A[x\[-UY ^ A[x][-2d} b ^ A[x][-rf] 4 ^ A[x 
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where the matrices of the differentials d{ and Ui are the same, % — 1, 2, 3, 4 (here 
again we set deg(x) = deg(y) = deg(^) = deg(w) = 1). We also consider the 
bi-graded Koszul complex on A[x\ associated to the sequence (x,y,z,w), and 
denote it (K m (x, y, z, w), v.) : 

-> A[x](-A) ^ A[x}(-3) 4 ^ A[x}(-2) 6 ^ A[x](-1) 4 ^ A[x\. 

The matrices of its differentials are obtained from the matrices of the differentials 
of by replacing f by x, fi by y, f 2 by z and f 3 by w. Note that since (x, y, z, w) 
is a regular sequence in A[x\, the complex K,(x,y, z,w) is acyclic. From both 
Koszul complexes (K,(f , fx, f 2 , fs), u.) and (K,(x, y, z, w),v,) we can build, as 
we did for curves, the approximation complex Z 9 . We define Zj := ker((i;) and 
Z,- t := Zi[id] (gu A[x\ for all i = 0, 1, 2, 3, 4 (where d : A — > 0), they are naturally 
bi-graded A [x] -modules. Since for all z = 1, 2, 3 we have Ui o t> i+1 + Vi o m+i = 0, 
we obtain the bi-graded complex : 

(Z.,v.) : - Z 4 (-4) ^ Z 3 (-3) ^ Z 2 (-2) ^ ^(-1) ^ Z = 

where, since we have supposed d > 1, Z4 — 0. 



Rema r k: In the la nguage of the moving surfaces method (see lSederberg and Chen 
|l995l |. Cox |200l| ). an element (gi, g 2 , g 3 , #4) G ^im(o) is nothing but a moving 



hyperplane of degree v following the surface S. 

To state the main result of this section, we need some notations. If p is an iso- 
lated base point defined by the ideal /, we denote by d p its geometric multiplicity 
(also called its degree); note that we have already denoted by e p its algebraic 
multiplicity. Recall that if M is a Z-graded .R-module, where R is a Z-graded 
ring, its initial degree is defined as indeg(M) = min{z/ G Z : M u 7^ 0}. 

Theorem 3.2: Suppose that the ideal I = (fo, fi, f 2 , f 3 ) C A is of codimension 
at least 2 and K is infinite. Let V := Proj(A/J) and denote by I-p the saturated 
ideal of I w.r.t. the maximal ideal m = (s,t,u) of A. Then, 

• The complex Z, is acyclic if and only ifV is locally generated by (at most) 
3 elements. 



Assume that V is locally generated by 3 elements, then for all integer 

v >v§:= 2(d — 1) — indeg(/p) 

the determinant of the graded complex of free K[x\-modules 

-> Z 3[u] (-3) ^ Z 2[v] (-2) ^ Z 1[u] (-l) ^ Z 0[u] = A [v] [x\ 

is an homogeneous element o/K[x] of degree d 2 — Yl P ep^p> an< ^ ^ s a multiple 
of independent ofv, where P is an implicit equation of S. It is exactly P 13 
if and only if I is locally a complete intersection. Moreover, for all v G Z, 
Z 3 ^ is always a free K[x]-module of rank max(( l '~2 +2 ) , 0); in particular 
Zsr u -\ = if and only if u < d — 1. 



L. Buse and M. Chardin: Implicitizing rational hypersufaces 



9 



Proof: See theorem 14.11 for the complete proof. Note that an argument similar 
to the one given in the proof of theorem 12.11 shows easily that Z% ~ d], 
and hence the last statement of the second point of the theorem. □ 

By standard prop erties of determinants of complexes (see e.g. appendix A in 
Gelfand etlrtl |l994j ) we deduce the 



Corollary 3.1: Suppose that I = (f Q , fi, f 2 , /s) C A is of codimension at 
least 2, K is infinite, and V = Proj(A/J) is locally generated by 3 elements. 
Then for all v > uq := 2(d — 1) — indeg(Jp) any non-zero minor of (maximal) 
size {y + 2)(z/ + l)/2 of the surjective matrix 

Zi\ v ]{-1) ^ A [u] [x] 
(gi, 92, 93, 9i) i-> xgx + yg 2 + zg 3 + wg 4 

is a non-zero multiple of P^ . Moreover, if V is locally a complete intersection, 
then the gcd of all these minors equals P". 

Before going through an algorithmic version of this theorem we make few 
remarks on the integer uq. First note that v§ depends geometrically on the ideal /, 
since indeg(J-p) is the smallest degree of a hypersurface in P 2 containing the closed 
subscheme defined by /. Let us observe how the initial degree of I-p behaves. If 
/ has no base points, that is I-p = A, then indeg(J'p) = 0. If there exists base 
points then this initial degree is always greater or equal to 1 since I-p is generated 
in degree at least 1, and is always bounded by d since the /j's are in l v . Also if 
/ is saturated then its initial degree is exactly d. We deduce that 

• Uq — 2d — 2 if I has no base points, 

• d — 2 < u < 2d — 3 if J has base points, 

• vq = d — 2 if I is saturated (note that in this case we know that Z3 i VQ i = 0, 
and hence det(Z,[ Uo j) is always obtained as a single determinant or as a 
quotient of two determinants). 

This shows in particular that the presence of base points simplify the complexity 
of the computation of the implicit sur face. Finally recall that the explicit bound 
(and not the theoretical one) given in Buse and Jouanolou 2002j | is only 2d — 2, 
if there exists base points or not. 



3.2. The algorithm 

We now develop the algorithm suggested by theorem 13.21 and then discuss some 
computational aspects. 

Algorithm (for implicitizing rational parametric surfaces with local almost 
complete intersection isolated base points, possibly empty): 
Input : Four homogeneous polynomials / 0; fi, fi, fz hi A of the same degree 
d > 1 such that the ideal (/o, /1, fi, f%) C A is locally generated by 3 elements 
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outside V(s, t, u) and at least of codimension 2. An integer v with default value 
v := 2d - 2. 
Output : Either : 

• a square matrix Ai such that det(Ai) is the determinant of (Z m )^, 

• two square matrices Ai and A2 such that ^[^2) is the determinant of 

• three square matrices A 1; A 2 and A 3 such that ^^wZ^f^ * s ^ ne deter- 
minant Of (2 u )t u y 

1. Compute the matrix F t of the first map of (J2J) : Af, — ^ A u+ d, and a kernel 
Ki of its transpose which has rank(Zi^) lines. 

2. Set m := iid^M±il_ Construct the matrix Zi defined by : 

Zi(i, j) = xKi(J,i) + yK 1 (j,i + m), zKi(j,i + 2m) + ioKi(j',i + 3m), 

with % = 1, . . . , m and j = 1, . . . , rank(^xM)- ^ JS ^ e m &trix of the map 
1) ^0 [1/] • 

3. If Zi is square then set Ax := Zi else 

(aj Compute a list L\ of integers indexing independent columns in Z t . L\ 
consists in m integers. Let Ai be the m x m submatrix of Z 1 obtained 
by removing columns not in L\ . 

(b) Compute the matrix F 2 of the second map of (J2J : A^ +d , and a 

kernel K 2 of its transpose which has rank(^ 2 [ iy j) lines, 
(cj Construct the matrix Z' 2 defined by, for all j — 1, ... , rank^^i), 

Z 2 (^j) = -yK 2 (:M) - zK 2 (j,i + m) - wK 2 (j,i + 3m),i = l,...,m, 
Z 2 (i, j) = xK 2 (j, i — m) — zK 2 (j, i + m) — wK 2 (j, i + 3m), i = m + 1, . . . , 2m, 
Z 2 (i,j) = xK 2 (j,i - m) + yK 2 (j,i) - wK 2 (j,i + 3m), i = 2m+l,.. .,3m, 
Z 2(hJ) = xK 2 (j,i) + yK 2 (j,i + m) + zK 2 (j, i + 2m), % = 3m + 1, . . . ,4m. 

(ci) Construct the rank(Zi^) x rank(2 2 [iy]) matrix Z 2 whose j th column 
Z2( # , j) J 's the solution of the linear system t Z 2 (»,j).K 1 = Z 2 («, j). It is 
the matrix of the map Z 2 \ u ](— 1) — ^> -Zim. 

(ej Define A' 2 to be the submatrix of Z 2 obtained by removing the lines 
indexed by L\. 

If A' 2 is square then set A 2 := A 2 eJse 

i. Compute a list L 2 of integers indexing independent columns in A' 2 . 
Define A 2 to be the square submatrix of A 2 obtained by removing 
columns not in L 2 . 

ii. Construct the matrix F 3 of the third map of (0) : ^4^ — ^ ^ +( i> an d 
the kernel K 3 of its transpose which has rank(Z^) lines. 
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Hi. Construct the matrix Z' 3 defined by, for all j = 1, . . . , rank(^3 







+ wK 3 (j,i + m),i = 1, ... 


,m, 




z' 3 (*,i) 


= -yK 3 (j, 


i — m) + u>K 3 (j, i + m),i -- 


= m + 1, . . 


.,2m, 






— 2m) + wK 3 (j, i + m),i = 


= 2m+l,. 


. , 3m, 




= -yK 3 (j, 


i — 2m) — zK 3 (j, i — m),i 


= 3m + 1, 


. . , 4m, 




= zK 3 (j,i 


— 3m) — zKs(j, i — m),i = 


4m + 1, . . 


. , 5m, 




= sK 3 (j',i 


- 3m) + yK 3 (j, i - 2m), i 


= 5m + 1, . 


. . , 6m. 



iv. Construct the rank^i^^]) x rank(^3[ jy ]) matrix Z 3 wiose j th column 
Z 3 («,j) is the solution of the linear system *Z 3 (»,j).K 2 = Z^*,^). 
It is the matrix of the map Z^uA— 1) ^2M- Define A 3 to be the 
square submatrix of Z 3 obtained by removing the lines indexed by 
L 2 . A 3 is of size 

endif 
endif. 

In the language of the moving surfaces method, the matrix Z x obtained at step 
2 gather all the moving hyperplanes of degree v following the surface S. Assume 
hereafter that V is locally a complete intersection. From a computational point 
of view, corollary 13.11 implies that this matrix can be taken as a representation 
of the surface S, replacing an expanded implicit equation (even if it is generally 
non-square). For instance, to test if a given point p = (xq : y : z : wq) E IP 3 is 
in the surface S, we just have to substitute x, y, z, w respectively by x , y , z , w 
in Zi and check its rank; p is on S if and only if the rank of Z 1 does not drop. Of 
course some numerical aspects have to be taken into account here, but the use 
of numerical linear algebra seems to be very promising in this direction. Note 
also that such matrices, whose computation is very fast, are much more compact 
representations of implicit equations compared to expanded polynomials which 
can have a lot of monomials. 

Even if the use of Zi seems to be, in the opinion of the authors, the best way 
to work quickly with implicit equations, the previous algorithm also returns an 
explicit description of P 13 , where P is as usual an implicit equation of S (always 
if V is locally a complete intersection). With the convention det(Aj) = 1 if Aj 
does not exists, for i = 2,3, P 13 can be computed as the quotient det ^ e 1 t ' ) ( .^ e 2 t | A3 ' > . 

Since det(Ai) is always a multiple of P 13 , we have det(Ai) = P^Q, where Q is an 
homogeneous polynomial in WL[x\. It can be useful to notice that this extraneous 
factor Q divides det(A 2 ); in fact det(A 2 ) = Qdet(A 3 ). 

3.3. Examples 

We illustrate our algorithm with four particular examples. It has been imple- 
mented in the software MAGMA which offers very powerful tools to deal with 
linear algebra, and appears to be very efficient. 
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3.3.1. An example without base points 
Consider the following example : 

fa = s% 
fl = t 2 u, 
fa = su 2 , 

fa = s 3+t 3 + U 3 . 

It is a parameterization without base points of a surface of degree 9. Applying 
our algorithm we find, in degree v = 2.3 — 2 = 4 a matrix Zi of size 15 x 24 which 
represents our surface. Continuing the algorithm we finally obtain a matrix A x 
of size 15 x 15, a matrix A 2 of size 9x9 and a matrix A 3 of size 3x3. Computing 
the quotient we obtain an expanded implicit equation : 

x 6 z 3 + 3x 5 y 2 z 2 + 3x A y 4 z + 3x A yz 4 + x 3 y 6 + 6x 3 y 3 z 3 + 3x 2 y 5 z 2 + 3x 2 y 2 z 5 — 
x 2 y 2 z 2 w 3 + 3xy 4 z 4 + y 3 z 6 . 

Trying our algorithm empirically with v = 3 we obtain a matrix Ai which is not 
square, it has only 9 columns (instead of the waited 10). The algorithm hence 
fails here, showing that in this case the degree bound 2d — 2 is the lowest possible. 

3.3.2. An example where the ideal of base points is saturated 



This example is taken from iSederberg and Chenl |1995| (see also iBuse et al 



[2002]), it is the parameterization of a cubic surface with 6 local complete inter- 
section base points : 

fo = s 2 t + 2t 3 + s 2 u + Astu + At 2 u + 3su 2 + 2tu 2 + 2m 3 , 

fi = -s 3 - 2st 2 - 2s 2 u - stu + su 2 - 2tu 2 + 2m 3 , 

f 2 = -s 3 - 2s 2 t - 3st 2 - 3s 2 u - 3stu + 2t 2 u - 2su 2 - 2tu 2 , 

f 3 = s 3 + s 2 t + t 3 + s 2 u + t 2 u — su 2 — tu 2 — u 3 . 

The ideal / = (/ , fi, fa, fa) is here saturated, so that we have u = d — 2 = 1, 
and we hence can apply the algorithm with v — 1. The matrix Zi is then square 
and the algorithm stops in step 2. It is given by : 

— z — w y + w 

x — 2y + z — 2w 2y — z 
—x — 2w y + 2w 

3.3.3. An example where the method of moving quadrics fails 

We here consider the example 3.2 in Buse et aD |2002j |. This ex ample was intro- 



duced to show how the method of moving quadrics (introduced in lSederberg and Chen 



[1995J) generalized in this paper to the presence of base points may fail. Consider 
the parameterization 

fa = su 2 , 

= t 2 (s + u), 
= st(s + u), 
= tu(s + u). 
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We can directly apply our algorithm in degree v = 2.3 — 2 = 4; we obtain a 
matrix Zi of size 15 x 30, matrices Ai and A2 of size 15 x 15, and a matrix A3 
of size 3x3. An expanded implicit equation is then xyz + xyw — zw 2 . 

If now we take into account that the ideal / = (fo, fx, fi, fs) has base points, 
we can apply the algorithm with v = 3. Even better, the saturated ideal of J, 
denoted I-p, is generated by t(s + u) and su 2 , showing that indeg(J-p) = 2 (since 
t(s + u) is of degree 2), and hence that we can apply the algorithm with v = 2. 
In this case Ai is 6 x 6, A 2 is 3 x 3 and the algorithm stops in step (e). 

Remark that to have the best bound possible for the integer v we have to 
compute the initial degree of I-p, which is not always obvious. However, it should 
be possible to determine some classes of parameterizations which satisfy certain 
geometric properties and for which we know in advance the initial degree of the 
saturation of /. 

3.3.4- An example with a fat base point 

With this example we would like to illustrate the theorem l3.2l when there is a base 
point minimally generated by three polynomials. Consider the parameterization 
given by 

f = s 3 - 6s 2 t - hst 2 - As 2 u + Astu - 3t 2 u, 
fx = -s 3 - 2s 2 t - st 2 - 5s 2 u - 3stu - Qt 2 u, 
f 2 = -As 3 - 2s 2 t + Ast 2 - 6t 3 + 6s 2 u - 6stu - 2t 2 u, 
f 3 = 2s 3 - 6s 2 t + 3st 2 - 6t 3 - 3s 2 u - Astu + 2t 2 u. 

The ideal I = {fo, fx, fi, h) defines exactly one fat base point p which is 
defined by (s,t) 2 . Therefore any implicit equation of our parameterized surface 
is of degree d 2 — e p = 9 — 4 = 5, and the degree of the determinant of the 
complex {Z,), , is of degree d 2 — d p = 9 — 3 = 6. Applying our algorithm with 
vq = 2(d — 1) — 2 = 2 we find that the matrix Z x is already square, of size 
6x6. Expanding its determinant one obtains a product of a degree 5 irreducible 
polynomial (an implicit equation of our surface) and a linear (since e p — d v = 1) 
irreducible polynomial. 

4. Implicitization of rational parametric hypersurfaces 

We know turn to the general problem of hypersurface implicitization, always 
using approximation complexes. In this section we prove new results to solve 
this problem under suitable conditions, and obtain the proof of theorem 13.21 as 
a particular case. 

We are going to consider hypersurfaces obtained as the closed image of maps 
from P n_1 to P ra , where n > 3. Consequently we set A := K[Xx, . . . , X n ], which 
will be the ring of the polynomials of the parameterizations we will consider. 
We also denote m := (Xx, ■ ■ ■ , X n ), and set — v := Hohia( — , A[— n\) and — * := 
Homgr A ( — ,A/m). Finally u;_ will denote the associated canonical module (see 
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Bruns and Herzog 1993]). Our first result is a technical lemma on the cycles of 



certain Koszul complexes. 

Lemma 4.1: Let f , . . . , f n be n + 1 polynomials of positive degrees do, ... , d n , I 
the ideal generated by them, K.(/; A) be the Koszul complex of the fi 's on A and 
denote by Zi and Hi the i-th cycles and the i-th homology modules of~K.(f; A), 
respectively. Denoting a := do + • • • + d n , if dim(A/ 1) < 1 then we have 

• Hi ^ for i = 0,1, Hi = for i > 2 and if 2 is zero if and only if 
dim(A/I) = 0. Ifdim(Afl) = l,H 2 ~ u A/I [n - a]. 

• If n > 3, then 



H m( Z P 



fori = 0,1 

H m( H i-p)*[ n ~ °1 fori = 2 

H*_ p [n - a] for 2 < i < n 

Z n- P l n ~ a \ f or i = n - 



Proof: The first point is classical, see e.g. iBruns and Herzo g 1993] 1.6.16 and 



1.2.4. For the second point we consider the truncated complexes: 

K> p : -> K n+1 -> ► K p+1 -> Z p -> 0. 

If p > 2 this complex is exact, and gives rise to a spectral sequence H^(K.> P ) 
which is at level 1 : 

■•• H° m (Z p ) 



■•• Hl-\Z P ) 

H m (K n+1 ) * H m (K p+1 ) * H m (Z p ), 

and modulo the identifications 

Hl{K.) ~ {K:r ~ (K'[-n]Y ~ (K n+1 _.[a - n])*, 
the last line becomes 

K*[n -a]^Kl[n-a]^-.-^ K* n _ p [n - a] - U£(Z P ). 
It follows 

• H l m (Z p ) = for i < p and for p + 2 < % < n, 
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. d> p „ p+1 :H*in-a]^HP m (Z p ), 
. d> p _ p :H*[n-a}^H^(Z p ), 

• dl^: H*[n- a] ^ H£\Z P ), 

where the upper index indicates from which spectral sequence the isomorphism 
comes, and the lower index indicates at which level of the spectral sequence the 
map is obtained. We also obtain a short exact sequence 

K_ p _ x [n - a] ^ K*_ p [n - a] -> H n m {Z p ) -> 

that gives the isomorphism H^(Z p ) ~ Z*_ p [n — a]. 

For p = 2, we have a spectral sequence if^(K> 2 ) =>- H^(H 2 ) which shows as 
in the previous cases that H^(Z 2 ) — Z*_ 2 [n — a], and that 

• H % m {Z 2 ) = for 4 < i < n, 

. d> 2 _ 2 :Hl[n-a]^HHZ 2 ), 
. d>l 3 :H* 2 [n-a]^H*(Z 2 ). 
It also provides an exact sequence 

- HHZ 2 ) ^ Hi{H 2 ) A - a] ^ H 2 m (Z 2 ) - 0, (3) 

where r is the transgression map of the spectral sequence. If H 2 ^ 0, we also 
look at the complex 

Kf 2 : -> Z 2 -> iT 2 -> ^ -> / -> 0, 

whose only homology is ifi and the corresponding spectral sequence H^(Kf 2 ) =>- 
H^(Hx). Noticing that = for i > 1, we get that h{(Z 2 ) = Hi(Z 2 ) = 

(this is also easily obtained by splitting K^ 2 into two short exact sequences, and 
taking cohomology). Together with (jSJ) and local duality, which gives H] n {H 2 ) ~ 
u>h ~ (A/ I sat )*[n — a] where 7 sat denotes the saturation of the ideal /, we get 
the asserted isomorphism for H^(Z 2 ). 

Note that we also get an isomorphism H^(H )*[n — a] ~ H^(Z 2 ) ~ H^(Hi), 
and an exact sequence: 

-> Hi{H x ) -> Hl{Z 2 ) -> f£(J) - 0. 

Finally for p = 1, the exact sequence 

K^ 1 : O^Z^Kx^I^O 

shows that H l - 2 (H ) ~ #™ V) - #m( z i) f o r * < ™, so that #°(Zi) = 
fl4(Z0 = 0, tf 2 ^) ~ H° m (H ) ~ - a], i^) ~ I^(tf ) ^ 

#*[n - a] and #4(Zi) = for 3 < i < n. 

The spectral sequence H 9 J(K> 1 ) (because only H^(H 2 ), 

and H^(Hi) may not be zero) gives H^(Zi) ~ Z^_ x [n — a], and this concludes 
the proof. □ 
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Remark: Note that the last spectral sequence of the proof also identifies H^(Zi) 
in another way by providing the exact sequence 

- Km T -U H{[n -a]^+ HliZ,) - 

which shows that if dim(A/7) = 1 

HjHliH^-UH^n-a]. 

Thus in this case Hi/ H^(Hi) has a symmetric free resolution. Also, the spectral 
sequences derived from the complexes K.f p provide exact sequences 

-> H ^) -> Hl{Z p ) -> /^(tf 2 ) - 

for p < n, and 

- - /7£ +1 (Z p ) - 77°(/7 ) - 

for p < n — 1, as well as an isomorphism H^ 2 (Z P ) ~ H^(H ) for p < n — 2. 

As we have already done in previous sections, we now consider approximation 
complexes. To do this, we introduce new variables T\, . . . , T n+ i which represent 
the homogeneous coordinates of the target of a given parameterization. Let 
/o, . . . , f n be n + 1 homogeneous polynomials in A of the same degree d > 1. 
Denoting by Zj the i th -cycles of the Koszul complex of the /j's on A, we set 
:= Zj[id] ®a ^.[Zj\ where [- -] stands for the degree shift in the JQ's and ( — ) 
for the one in the Tj's. It appears that the differentials t>, of the Koszul complex 
K,(Ti, . . . ,T n+ i, A[T\) induce maps between the Z^s, and hence we can define 
the approximation complex (note that Z n+ i = 0) 

(Z„ v.) : -> Z n (-n) ^ . . . ^ Z 2 (-2) ^ ^(-1) ^ Z = A[T]. 

It is a naturally a bi-graded complex, and it is easy to check that H (Z,) ~ 
Sym A (7). We now give some acyclicity criterions in case / = (/ , . . . , f n ) define 
isolated points in Proj(A). Recall that a sequence xi, . . . ,x n of elements in a 
ring R is said to be a proper sequence if 

Xi + \Hj(xi, . . . , xf, A) — for i = 0, . . . , n — 1 and j > 0, 

where the Jf/s denote the homology groups of the associated Koszul complex. 

Lemma 4.2: Suppose that I = (/ , . . . , / n ) is of codimension at least n — 1, K 
is infinite, and let V := Proj(A/7). Then the following are equivalent : 

(1) Z, is acyclic, 

(2) Z, is acyclic outside V(m), 

(3) I is generated by a proper sequence, 
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(4) V is locally defined by a proper sequence, 

(5) V is locally defined by n equations. 



Proof: Bv iHerzog etaD [l983j , (1)^(3) and (2)^(4). Moreover (1)=^(2) and 



(3)=j>(4). We will show that (4)^(5)^(3). 

Assume (4). For each p G V there exists a non empty open set fl p in K( n+1 ) 
such that if (ay) G fi p and g>j := Ylj a ijfji the g^'s form a proper sequence 
locally at p. Taking (ay) G : = n pG pJl p , this gives g±, . . . ,g n +i so that the g^'s 
form a proper sequence outside V(m). We may also assume, by shrinking Q if 
necessary, that gi, . . . , <? n -i is a regular sequence, so that gi,...,g n is clearly 
a proper sequence. Now g\,...,g n+ i is proper if and only if g n+ i annihilates 
Hi(gi, . . . ,g n ; A) ~ oja/j (locally outside V(m), a priori), where J denotes the 
saturated ideal of (g±, . . . ,g n ) w.r.t. m. But Ami^^A/j) = J, so that g n+ \ G J, 
which proves (5). 

Now assume (5). One considers (in the same spirit as above) g%, . . . , g n+ i so 
that gi, . . . , g n -\ is a regular sequence and g%, . . . , g n defines V . Then g n+ i G I-p, 
and therefore annihilates ua/i v , so that the g^s forms a proper sequence. □ 

We are now able to state our main result on the hypersurface implicitization 
problem. 

Theorem 4.1: Let n > 3. Let I be the ideal (f , . . . , f n ), V := Proj(v4/7) and 
I-p the saturation of I w.r.t. m. Assume that IK is infinite and that dim'P < ; 
that is V define a finite number of base points in Proj(v4), possibly empty. Then 
we have 

• Z, is acyclic if and only if V is locally defined by n equations. 

• Let v > v Q := (n — l)(d — 1) — indeg(Jp). Then H^(Sym A (I))[u] = 0. 
Moreover (Z m )[ v ] ®k[t] K(T) is acyclic if and only if Z, is acyclic. 

• Let v > u and assume that V is locally defined by n equations. Then 
D := det ((£,)[„]) is a non zero homogeneous element o/K[T], independent 
of v (modulo K. x ), of degree d 71 ^ 1 — ^2 p&v d p . Denoting by 7i the closed 
image of the rational map : P n_1 — > P n ; D = H dcg ^G where H is an 
implicit equation ofTL. Moreover G G K x if and only ifV is locally of linear 
type, and if and only ifV is locally a complete intersection. 

Before giving the proof of the theorem we recall that, by defi nition, V is said to 



be loc ally of linear type if Proj(Sym A (/)) = Proj(ReesA(/))- Bv lBuse and Jouanolou 
|2002j | theorem 2.5, deg(H deg ^) = d n ~ l — ^ p( z V e p , and consequently we always 
have deg(G) = J2 P ev( e p ~ ^p)- ^lso e p — ^p with equality if and only if V is 
locally a complete intersection at p G V. 

Proof: The first point follows from lemma 14.21 For the second point we con- 
sider the two spectral sequences associated to the double complex H^(Z 9 ), both 
abouting to the hypercohomology of Z,. One of them abouts at level two with: 
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H^H^Z.)) for p = 0, 1 and q > 
2 >E P = 00 >E P ={ #&(Sym A (I)) for g = 

else. 

The other one gives at level one: 

i"E P q = R p m {Z q )[qd} ® A A[T]{-q). 

By lemma 14.11 H^(Z q ) = for p < q and for p < 2. This shows that 
H^(H q (Z,)) = for q > except possibly for p = q = 1, by comparing the 
two spectral sequences. Therefore, H q (Z.) = for g > 1 and If (.£,)) = 0. 

If p > 2 

~ H*[n - (n + 1 - p)d\ ® A A[T](-p) 
so that {i"E p p ) v = if v > (n - 2)d - n. Also, 

~ (/ P //)*[n - (n - ®a A[T](-2), 

so that (i'EDv = for z/ > (n — l)c? — n — indeg(/p/J). 

Therefore, if z/ > z/ , (i/'-Epj, = for p < q. Note also that we have the 
equalities 

mm{d, indeg(/p//)} = min{<i, indeg(Jp)} = indeg(Jp). 

By comparing with the other spectral sequence, we have H^Hi^Z,))^ = and 
H° m (Sym A (I)) v = 0. 

As if^(Sym A (/)) i/ = for v > b> we have 

Ann K [T](Sym A (/) I/ ) = Ann K [T](Sym A (/) i , ) 

for any v > u (see for instance iBuse and Jouanoloul [2002] , the proof of 5.1), 
so that this module is torsion if and only if V is locally defined by at most n 
equations (because AnnK[r](Sym A (/)y) is torsion for v ^> if and only if / is 
defined by < n + 1 equations outside VCm); one may a lso use the study of the 



minimal primes of Sym A (J) in lHuneke and Rossil |l986j |). Also in the case where 



V is locally defined by at most n equations, the divisor associated to this module 
is independant of v > Uq- This finishes the proof of point two. 

We come to the third point. Note that we have just shown that D is indepen- 
dant of v (up to an element of K x ). To compute the degree of D, we may then 
take v 3> 0. 

The matrices of the maps of Z m have entries which are linear forms in the Tj's, 
so that the determinant of {Z,), , is a form in the Tj's of degree 



<5:=^(-l) i+1 ^dim K (^ [l/+ld] ). 



i=l 



L. Buse and M. Chardin: Implicitizing rational hypersufaces 



19 



We have canonical exact sequences, i = 0, . . . , n, 

-> -> ^ -»■ 

and 

-> 5, -> -> -> 

which shows that S can be expressed in terms of the Hilbert polynomials of the 
Ki's and the H^s. The contribution of the K^s only depends on n and d, and 
the one of the H^s only comes from Z\ and Z 2 and is : 



2d 



as degi^i = 2deg"P (one may use for example that (H ) u — (Hi) u + (H 2 ) u = 
for v ^> and that H 2 — ojr/i v up to a degree shift), this contribution is equal 
to — degP for v ^> 0. The contribution of the K[s is cP -1 (for v ^> 0) in the 
case where the H^s are for v ^> 0, and we get 5 = d n ~ x — degP. 

Let X := Proj(Rees 7 (A)) C K := Proj(Sym A (J)) C P 71-1 x P n . The scheme 
X is the closure of the graph of (f) in P ra_1 x P n . Therefore, the closed image 
of <p is p 2 (X), and p 2 (X) C p 2 (y). As D ^ 0, p 2 (F) ^ P n . If V is locally of 
linear type, then X = Y, so that G G K x which implies that d p = e p for any 
p G V and therefore "P is locally a complete intersection. If V is not locally a 
complete intersection at p G V, d p ^ e p so that deg G > and 1^7. Note that 
these last f acts also follows from the study of the minimal primes associated to 



Sym A (I) in lHuneke and Rossi 



1986. □ 



This theorem yields easily an algorithm to implicitize parameterized hyper- 
surfaces under suitable assumptions, using only linear algebra routines. The case 
n = 3 have been completely explicited in section El 



References 

Winfried Bruns and Jurgen Herzog. Cohen-Macaulay rings. Cambridge studies 
in advanced mathematics, 39, 1993. 

✓ 

Laurent Buse. Etude du resultant sur une variete algebrique. PhD thesis, Uni- 
versity of Nice, 2001. 

Laurent Buse, David Cox, and Carlos D'Andrea. Implicitization of surfaces in 
P 3 in the presence of base points. Preprint math. AG/0205251, 2002. 

Laurent Buse and Jean-Pierre Jouanolou. On the closed image of a rational map 
and the implicitization problem. Preprint math. AG/0210096, 2002. 

David A. Cox. Equations of parametric curves and surfaces via syzygies. Con- 
temporary Mathematics, 286:1-20, 2001. 



L. Buse and M. Chardin: Implicitizing rational hypersufaces 



20 



David A. Cox, Ronald Goldman, and Ming Zhang. On the validity of implic- 
itization by moving quadrics for rationnal surfaces with no base points. J. 
Symbolic Computation, 29, 2000. 

David A. Cox, Thomas W. Sederberg, and Falai Chen. The moving line ideal 
basis of planar rational curves. Comp. Aid. Geom. Des., 15, 1998. 

Carlos D' Andrea. Resultants and moving surfaces. J. of Symbolic Computation, 
31, 2001. 

I.M. Gelfand, M.M. Kapranov, and A.V. Zelevinsky. Discriminants, Resultants 
and Multidimensional Determinants. Birkhauser, Boston-Basel-Berlin, 1994. 

J. Herzog, A. Simis, and W.V. Vasconcelos. Koszul homology and blowing-up 
rings. Lecture note in Pure and Applied Math., 84:79-169, 1983. 

Christoph M. Hoffmann. Geometric and solid modeling : an introduction. Mor- 
gan Kaufmann publishers, Inc., 1989. 

C. Huneke and M.E. Rossi. The dimension and components of symmetric alge- 
bras. J. of Algebra, 98:200-210, 1986. 

J.-P. Jouanolou. Resultant anisotrope: Complements et applications. The elec- 
tronic journal of combinatorics, 3(2), 1996. 

T.W. Sederberg and F. Chen. Implicitization using moving curves and surfaces. 
Proceedings of SIGGRAPH, pages 301-308, 1995. 

A. Simis and W.V. Vasconcelos. The syzygies of the conormal module. American 
J. Math., 103:203-224, 1981. 

W.V. Vasconcelos. Arithmetic of Blowup Algebras, volume 195 of London Math- 
ematical Society Lecture Note Series. Cambridge University Press, 1994. 



